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Abstract 



o 

■ Context. CMB experiments aiming at a precise measurement of the CMB polarization, such as the Planck satellite, 

need a strong polarized absolute calibrator on the sky to accurately set the detectors polarization angle and the cross- 
polarization leakage. As the most intense polarized source in the microwave sky at angular scales of few arcminutes, 
the Crab nebula will be used for this purpose. 
Aims. Our goal was to measure the Crab nebula polarization characteristics at 90 GHz with unprecedented precision. 
Methods. The observations were carried out with the IRAM 30m telescope employing the correlation polarimeter XPOL 
and using two orthogonally polarized receivers. 
Results. We processed the Stokes /, Q, and U maps from our observations in order to compute the polarization angle and 
linear polarization fraction. The first is almost constant in the region of maximum emission in polarization with a mean 
value of asky = 152.1 ± 0.3° in equatorial coordinates, and the second is found to reach a maximum of II = 30% for the 
most polarized pixels. We find that a CMB experiment having a 5 arcmin circular beam will see a mean polarization 
angle of asky = 149.9 =b 0.2° and a mean polarization fraction of II = 8.8 =b 0.2%. 



g : Key words. ISM: supernova remnant, Polarization, (Cosmology:) cosmic background radiation 
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1. Introduction Crab nebula as an absolute polarization calibrator for cos- 

mic microwave background (CMB) experiments, such as 

The Crab nebula (Tan A, Ml or NGC1952, at coordinates the planck EgA mission . The results for circular polar _ 

a = 5 34 m 32 s and S = 22°0 52 , J2000) is a supernova ization are irrelevant for current CMB p o larimetry exper- 

remnant that emits a highly polarized signal due to both iments a nd wil l be published elsewhere (Wiesem ever^tld] 

the synchrotron emission of the centra l pulsar and i ts inter- |j n preparation! ) 

action with the surrounding gas (see eg Hester 2008). In this 

paper, we present the measurements of intensity and polar- The last Thomson scattering of cosmological photons 
ization of the Crab nebula at 90 GHz (3.3 mm) performed by free electrons at the time of decoupling; between mat- 
at the IRAM 30 m telescope using the XPOL instrument ter a nd radiation at z dec = 1088.2 ± 1.1 (|Dunklev et al.l 
(|Thum et al.ll2008[ ). The strength of this instrument is the l2QQ9h have led to a linearly polarized CMB radiation at 
precision in the angle of polarization with respect to the a level of a few percent. The linear polarization pattern 
sky reference frame, which highly improves previous mea- on the sky can be decomposed into a scalar and a pseudo- 
surements. This precision is needed in particular to use the scalar field, respectively called E- and 5-modes, which have 
the advantage of being independent of the re ference frame 



tne ad vantage oi being independent ol tne r eference irame 
(jZaldarriaga fc Seliaklll997HZaldarriagal200l[ ). This decom- 
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fluctuations (scalar perturbations) can only produce E po- 
larization, while B can only be produced by lensing ef- 
fects, exotic pheno mena such as cosmological birefringence 
([Carroll et al.lll990l ) and by the stochastic background of 
gravity waves arising from the inflation epoch (tensor per- 
turbations). The detection of the latter would improve the 
constraints on the inflationary model parameters in an in- 
valuable manner and constitutes one of the most ambitious 
goals of observational cosmology. 

The first measurement of the CMB E polarization was 
made by the 30 G Hz radio interferometer, DASI, in 2002 
([Kovac et al.l [2002). Since then, in addition to a furthe r 
measurement by the DASI experiment ([Leitch et al.ll2005[ ), 
E modes measuremen ts have been made with the CB I 
(jReadhead et al.l |2004ft. CAPMAP ([BischofT et al.l 1 2008). 
BOO MERANG (IMontrov et al 1120061 ). Wmap (iPage et al 



20071) . MAXIPQ LdWu et all I2007D . QUaD ( Ad e et al 



2008D and Bicep (Chi ang et al.ll2QQ9h . With these measure 



ments, the Cf E angular power spectrum has now been de- 
tected over a wide range of angular scales (30 < £ < 1000). 
Additionally, the Wmap satellite has given a very pre- 
cise measurement of the temperature- polarization cor- 
relation, the Cj E angular cros s power spectrum , at large 
scales up to multipoles £ ~ 500 (iNolta et al.ll2009D . To date, 
no detection of non-zero Cf B angular power spectrum has 
been reported and only upper limits on the effect of gravity 
waves from other spectra can be inferred. The best indirect 
limit to date on the tensor to scalar ratio r, given by Wmap 
data and ot her cosmo logical obser vables together is r < 0.2 
(95% C.L.) (jKomatsu et al.ll2009[ ). 

Upcoming experiments target the detection of lower val- 
ues of the tensor-to-scalar ratio. The low signal of the B- 
modes compared to the temperature and £?-modes signals 
will require a high control of systematics and good knowl- 
edge of the instrument. In particular, the knowledge of its 
polarization characteristics will require a highly and well 
known polarized source. 

The Crab nebula represents the most suitable candidate 
as absolute calibrator for the polarization angle and the 
linear polarization degree, given that it is the most intense 
polarized source in the microwave sky at angular sca les of 
few arcminutes ([Flett fc Henderson|[l979l Il98ll Il983l ). 

The paper is organized as follows: in §[2] we describe the 
required calibration accuracy for a CMB experiment aiming 
at the detection of the 5-modes; in §[3] we briefly describe 
the XPOL instrument; § |4] reports the observations, data 
reduction and systematic error analysis; in §[5] we give our 
results; we discuss the extension of our measurements to 
other frequencies in § [6] and summarize our results in § [7J 

2. Required calibration accuracy 

The expected primordial B-mode level at large angular 
scale, for r = 0.01, is typically 1% of the I?- mode level. More 
precisely, for £ < 20, the E-mode is on t he order of 0.07 /iK 2 
(from the WMAP best-fit model, see Nolta et al.l 120091 ). 
while the B-mode is expected at 10 _4 /iK 2 . For compari- 
son, the temperature anisotropics are at the level 1000 /iK 2 
at large scale. The Planck High Frequency Instrument, for 
example, has the potential to detect or constrain the tensor- 
to-scalar ratio down to r ~ 0.03 (Efstathiou & Gratton 
2009). However, such a low level will be reachable only if the 
systematics are well controlled, otherwise, given the hierar- 
chy of signal levels, T ^> E ^> £?, we may expect leakage 



from temperature to polarization and a leakage of .E-mode 
to B-mode. 

Various instrument parameters, if not precisely known, 
will induce systematic effects. In the case of the Planck HFI 
polarization sensitive bolometers, the detectors are sensi- 
tive to the power carried by the electric field aligned with 
the detector's grid. The signal s measured by a detector 
can be written ([Rosset et al.ll2010[ ) 



9 



I 



I + p(Qcos2(a + a d ) 



U sin 2(a + ad)) 



A(n)dn, 



(i) 



where /, Q and U are the Stokes parameters of the sky sig- 
nal (intensity and linear polarization in the sky reference 
frame), g is the detector gain, A is the instrument beam, p is 
the detector polarization efficiency, a is the angle between 
the telescope reference frame and the local meridian in the 
sky, ad is the angle at which the polarimeter is maximally 
sensitive to polarized radiation, in the telescope reference 
frame (here, we have made the simplification that the beam 
for polarization is identical to the intensity beam, only mod- 
ulated by the sine and cosine of the detector orientation) 
and n is a direction on the sky. In this simplified model, 
g and A are measurable with unpolarized light (cosmologi- 
cal and orbital dipole, planets). Only the polarization spe- 
cific parameters, a and p, require a known polarized signal 
to be calibrated. Measuring with precision these proper- 
ties in labor atory or with an artincal source can be d iffi- 
cult (see ea iRosset et~aD 120101 : fTakahashi etTaD I2009D . so 
measuring them with a known astrophysical source allows 
cross-checking or even a better calibration. In any case, the 
best limit on the precision of the polarization efficiency and 
the orientation of detectors will be set by the precision on 
the polarization fraction and polarization orientation of the 
source. 

In the case of a unique detector scanning the s ky, it 
has been shown (|Q'Dea et gil.ll2QQ7c IRosset et al.ll2010[ ) that 
the error on the B-mode power spectrum due to an error 
Aa in the polarization orientation and an error Ap on the 
polarization efficiency is given by: 



ACf B = 2ApCf B + 4Aa 2 Cf E 



(2) 



For example, an experiment targeting a tensor-to-scalar ra- 
tio measurement of r = 0.01, using the large angular scale 
B-mode (at £ < 20), will require the polarization orien- 
tation and the polarization efficiency to be known with a 
precision of: 



Aa < 0.9°, Ap < 0.05, 



(3) 



in order to have a leakage from E to 5-modes lower than 
one tenth of the expected Cf B . 

These tolerances set the precision needed for the cali- 
brator, and define our goal for the Crab nebula observations 
presented in this paper. 

3. Instrumental setup 

3.1. General 

For the observation of the Crab nebula, we used the 
IRAM 30 m telescope at 89.189 GHz, the frequency of 
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the HCO + (1-0) transition, where the angular resolution 
(the full width at half power of the near-Gaussian beam) 
is 27". We used the following IRAM receivers: A100 cryo- 
stat (vertical linear polarization with respect to Nasmyth 
reference frame) and B100 cryostat (horizontal linear polar- 
ization). The mean system temperatures were 87 K (A100) 
and 113 K (B100) while mean receiver temperatures (single- 
sideband) were 67 K (A100) and 59 K (B100). The ef- 
fective cold loads used for temperature calibration had a 
mean temperature of 87 K (A100) and 98 K (B100). We 
used the VESPA backends in XPOL mode (see §|33, with 
a 500 MHz bandwith (200 channels spaced by 2.5 MHz). 
Finally, the conversion factor from antenna temperature 
outside atmosphere to flux density was 1 K = 6.0 Jy and 
the mean zenith opacity was r z =0.06 at our observing fre- 
quency. 

3.2. XPOL 

The signals from the two orthogonal linearly polarized 
heterodyne receivers were detected in auto- and cross- 
correlation, from which the four Stokes parameter s were 
derived as described in detail by iThum et al.l ((2008). This 
procedure, designated XPOL, employs a precise calibration 
of the phase between the two receivers, resulting in a phase 
error < 1° per channel of 2.5 MHz, across the 500 MHz 
bandwidth. 

Particular care was taken in the absolute calibration 
of the polarization angle r. This angle was re-measured, 
for the purpose of this experiment in the Nasmyth cabin, 
with respect to its horizontal axis, where the receivers are 
stationary, as 

r-i«Mg), (4) 

where Q and U refer to the quantities measured in the 
Nasmyth reference frame. The precision with which the re- 
ceivers, labelled H and V, actually measure horizontally 
and vertically polarized power is dete rmined by t he or i- 
entation of the grid (G3 in Fig. 1 of lThum et ail (|2008h ) 
which splits the incoming beam into a transmitted frac- 
tion which is horizontally polarized, and a reflected frac- 
tion which is vertically polarized. The horizontality of the 
incoming beam, and thus also that of the H-fraction rests 
on the correct alignment of the telescope Nasmyth mirror 
M3, the subreflector, and the validity of the implemented 
homology corrections and the pointing model. Their com- 
bined error is less than 1 arcmin. The verticality of the 
V-fraction additionally depends on the orientation of G3 
whose normal should make an angle of 45° with respect to 
the local vertical. The mounting frame of G3 has a typical 
mechanical precision of ^ 1 mm over the length of 25 cm 
of the grid support frame. This corresponds to an angular 
precision of ~ 0.3° which was confirmed by a direct mea- 
surement with an electronic angle gauge in contact with the 
grid's frame. 

Equation (jH shows that the precision of the r measure- 
ment further depends on the correct relative calibration of 
Stokes U and Q. Since XPOL measures Q as the power dif- 
ference between the H and V receivers, but derives Stokes 
U from their correlation, we need to make a precise mea- 
surement of the inevitable correlation losses. This is done 
as a by-product of the frequent phase calibration measure- 



ments where a wire grid of precisely known orientation is 
observed against hot and cold loads. Due to the high signal- 
to-noise of these calibrations, the correlation losses can be 
determined with an r.m.s. error of ~ 1%, resulting in an 
r.m.s. error of r of ~ 0.3°. 

The transfer of r to the polarization angle asky as mea- 
sured from north to east in the equatorial system is made 
through 

«Sky = 90 o +r-(e-7?), (5) 

where the elevation e and the parallactic rotation angle rj do 
not introduce any additional error. We therefore estimate 
that our polarization angles have an absolute precision on 
the order of 0.5° or slightly better. 

4. Observations and data reduction 

4.1. Mapping strategy 

The Stokes maps of the Crab nebula were taken in the 
on-the-fly mode (hereafter OTF), allowing us to map a 
large area on the sky within a reasonable time interval. 
The telescope was scanning the source while XPOL dumped 
twice per second the total power from the orthogonally po- 
larized receivers, and the complex cross-correlation. The 
Stokes maps were subsequently derived from these quanti- 
ties. Throughout the paper, Q and U will refer to the quan- 
tities measured in the Nasmyth system and then rotated in 
equatorial coordinates according to equation The Crab 
nebula was scanned along the projections of latitude and 
longitude onto the local tangential plane (hereafter respec- 
tively a and 6 scanning directions). Such a scanning strat- 
egy is indispensable if artifacts due to the variations of the 
atmospheric emission, leading to stripes in the maps, are 
to be removed (see section 4.2). The maps scanned along 
the a direction are 8.8' x 6.4' large, those scanned along 
the 5 direction 10.2' x 5.5'. The differences in size are due 
to slightly different system overheads coming from the ac- 
celeration and deceleration of the antenna mount drive be- 
fore and after scanning an OTF stripe. The scanning speed 
was 16"/sec, leading to a smearing of the 27" wide beam 
(FWHM) by 8". A pair of maps scanned in the orthogonal 
directions took 22 minutes, including regular temperature 
and phase calibrations. All the maps were projected on the 
same 51x41 pixels frame with square pixels of size 13.7". 

4.2. Raw data reduction 

The conversion from backend count rates to antenna tem- 
peratures follows the standard method applied at the IRAM 
30 m observatory, involving total power measurements of 
the sky, a hot load, and a cold load. These calibration mea- 
surements are combined with the phase calibration mea- 
surements, namely an on-off between a fully linearly polar- 
ized signal at cold load temperature, and the unpolarized 
hot load signal at ambient temperature (§3). The calibra- 
tions were de rived and applied to the data with the MIRA 
software (see Wie semeyerl T2009) . 

Since the continuum emission from the Supernova rem- 
nant has a well-defined outer boundary, there is no need 
to observe a separate emission-free reference position. The 
atmosphere emission was rather removed by subtracting, 
from each OTF subscan observed along a given direction, 
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a linear baseline defined by the Crab nebula emission-free 
map boundary. Subtracting only a zero order spatial base- 
line would leave "stripes" in the maps along the scanning 
direction, which co uld be removed in Fourier space with the 
PLAIT algorithm (H merson fc Graevdfl988h . Both meth- 
ods were shown to yield equivalent results. In fact, they are 
both limited by the uncertainty of the absolute reference 
flux if the source boundary is not precisely known. In the 
image plane approach, the source and sky emission cannot 
be clearly separated from each other, in the Fourier plane 
technique case, the PLAIT algorithm does not work close 
to zero spatial frequency (whose Fourier component is the 
total flux density contained in the map). 

Since for XPOL the signal is downconverted into eight 
adjacent spectral basebands of 62.5 MHz width each, to- 
tal power variations between the sky/hot/cold scan and 
the OTF record result in calibration mismatches among 
the basebands that introduce platforms in the autocorre- 
lation spectra. In an attempt to correct for this undesired 
effect, one may naively think of using the data from a fil- 
terbank, connected in parallel and with a monolithic spec- 
tral band, such that they are not affected by platforming. 
However, in the subsequent data processing step (rotation 
from the Nasmyth to equatorial reference frame, for getting 
a stationary astronomical polarization), the spectra will be 
linearly combined with the real part of XPOL's cross cor- 
relation (which is, after phase calibration, Stokes U in the 
telescope's Nasmyth reference frame). The small calibration 
differences observed between XPOL and the filterbank are 
believed to be caused by non-linearities when the backends 
operate close to the saturation limit. Mixing these data 
would thus lead to an inconsistent calibration and hence 
an error in the polarization angle map of the Crab nebula. 
Therefore, another strategy was used: 

— For each autocorrelation spectrum (i.e. each record of 
an OTF map), the median across the whole band, with 
concatenated basebands, has been calculated. 

— For each of the eight basebands of a given spectrum, the 
mean continuum flux and noise has been determined. 

— When the absolute difference between the mean base- 
band signal and the median signal is greater than twice 
the baseband noise, the offset of the baseband signal 
has been corrected by the difference between the me- 
dian and the baseband mean. 

The result is a spectrum that looks "flat" . The remain- 
ing calibration mismatches lead to stripes in the maps, as 
do fluctuations of atmospheric origin, and gain fluctuations 
owing to receiver instabilities. As a matter of fact, stripes 
are also present in filterbank spectra, which do not suffer 
from platforming. 

It should be noted here that the cross-correlation spec- 
tra (Stokes U in the Nasmyth reference frame, and Stokes 
V) do not suffer from platforming either. This derives 
from the fact that the clipping voltage of a sampler is ad- 
justed with respect to the signal's noise level, which may 
change when a given interval of the time series is mea- 
sured. The zero time lag channel measures the total power 
(Parseval's theorem), which is positive definite. After the 
FFT, this shows as an offset to the spectral baseline. A 
cross-correlator analyses the signals from two independent 
samplers. Variations in the noise power of the signals com- 
ing from the samplers are largely uncorrelated, and there- 



fore cancel out in the spectral cross-power (which is not a 
positive definite quantity anymore). 

The data were gridded at Nyquist sampling. Because of 
the slightly different values of the start or end scanning co- 
ordinate between different subscans, the sampling is irreg- 
ular. Therefore the data had to be resampled to a regular 
grid. This has been done by applying a Gaussian convo- 
lution kernel with a FWHM of one-third the telescope's 
half-power beam width. The kernel is truncated at three 
FWHM. 

4.3. Description of the data sets 

We performed two observation campaigns of the Crab 
Nebula, respectively from September 5 to 10, 2007, and 
from January 9 to 12, 2009. 

— First campaign: Observations were made during 6 inter- 
vals, under varying astmospheric opacity and stability 
conditions. Most of the maps were affected by a sub- 
stantial atmospheric contamination, showing as strong 
linear stripes in the scanning direction. We obtained 66 
/, Q, U and V maps resulting in 8' x 6.4' maps after 
coadding a-scanned and (5-scanned maps. These maps 
were centered on the reference point we called crabx- 
pol having equatorial coordinates a = 5 h 34 m 31.5 s and 
5 = 22°00'27.7" (J2000). The centering and the size 
of these maps made the destriping delicate because the 
extension of the source emission was larger than the ob- 
served area and thus, strong mix-up with atmospheric 
contamination arose on maps' edges leading to uncon- 
trolled systematics. 

— Second campaign: Observations were made during 4 
nights during which we experienced stable atmospheric 
conditions. As a consequence of the first observation 
campaign, we increased the coaddition of the a-scanned 
and ^-scanned maps to a size of 10' x 6.7' in order to be 
able to characterize the atmospheric stripes outside the 
source. 11 maps of /, Q, U and V were observed centered 
on the same reference point as in the first campaign, 
namely crabxpol, and then we changed to a new refer- 
ence point for center, called crabcenter and having equa- 
torial coordinates a = 5 h 34 m 31.5 s and S = 22°01 / 17.7 // 
(J2000), on which the next 28 maps were observed. This 
choice was made to have a better centering of the source 
extent in the north-south direction, allowing better de- 
striping. The 11 first maps centered on crabxpol have a 
significantly higher noise level than the 28 centered on 
crabcenter due to worse weather condition at the begin- 
ning of the campaign. 

These data sets thus gather 105 individual J, Q, U and 
V maps, divided in three subsets: first campaign maps cen- 
tered on crabxpol (hereafter 1CX), second campaign maps 
centered on crabxpol (hereafter 2CX) and second campaign 
maps centered on crabcenter (hereafter 2CC). Because in 
the two first subsets the source extends over an area larger 
than the mapped one, the destriping is delicate and thus 
may lead to uncontrolled systematics. For each subsets 
and for combinations of different subsets we computed, for 
each Stokes parameter, mean and dispersion for each pixel 
among the whole data set. As the pixels are small with re- 
spect to the beam, the computed dispersion does not suffer 
from bias. Dispersion maps, in the case of the first and sec- 
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(J (I) Crab, all 39 maps <j(I) Crab, crabcenter maps 




4 2 -2 -4 4 2 -2 -4 

Aa ( arcmin ) Aa( arcmin ) 



Figure 1. Standard deviation maps for the Crab nebula intensity a (I) (top row) and polarized intensity <j(P) (bottom 
row), in antenna temperature (K). a (I) and cr(P) maps were computed using the 2CC and 2CX (see text) 39 maps 
(left column) and using the 28 2CC maps only (right panel). / and respectively P contours computed with the 28 2CC 
maps are overploted. We can see that when maps scanned along the a direction are added to the data set the standard 
deviation is higher and shows some features indicating an uncontrolled systematic effect. 



ond subsets, present specific high-dispersion features aris- 
ing from the confusion between stripes and source emission 
on the borders during the destriping procedure. This can 
be seen on Figure [T] where the combination of subsets 2CX 
and 2CC, for the intensity, shows a higher standard devia- 
tion feature than the 2CC alone. For consistency, we also 
checked the behavior of the polarized intensity standard 
deviation which slightly decreases when using only 2CC 
maps. 

Thus, using only the 28 maps from 2CC provides the 
best results and any addition of 2CX maps, only 2CX maps 
scanned in a direction (direction in which the borders of the 
maps are off-source) and 1 CX lead to worse results in terms 
of noise homogeneity and of signal-to-noise ratio, both in 
intensity and polarization. Consequently, our analysis will 
use only the 28 maps of the second campaign centered on 
the crabcenter reference point. 

4.4. Systematic effects analysis 

In order to check for residual systematic effects in the maps, 
we performed jack-knife tests dividing the 28 maps in sev- 
eral subsets using three criteria, i.e.: randomly, between 
maps scanned in a and S direction and between the three 
observation nights during which these 28 maps were taken. 



Data subset 


Pol. angle asky 


Pol. fraction II 


Random 1 
Random 2 


152.19 ±0.36° 
152.07 ±0.30° 


23.77 ± 0.52% 
23.91 ± 0.54% 


a scans 
6 scans 


152.02 ±0.32° 
152.26 ±0.35° 


23.38 ± 0.45% 
24.41 ± 0.53% 


Night 1 
Night 2 
Night 3 


152.47 ± 0.37° 
151.89 ±0.31° 
152.10 ±0.30° 


23.96 ± 0.55% 
23.93 ± 0.56% 
23.56 ± 0.49% 


XoG{-30°,45°} 
Xoe{45°,120°} 


152.17 ±0.37° 
152.09 ±0.30° 


23.60 ± 0.52% 
24.12 ±0.55% 


Whole data set 


152.13 ±0.34° 


23.84 ± 0.54% 



Table 1. Polarization angle asky and polarized fraction II 
computed for different divisions of the data set. 



We also divided the data in two with respect to the posi- 
tion angle of the Nasmyth reference frame projected onto 
the sk y, yp, in order to ch eck for sidelobe polarization ef- 
fects ( Forbrich et al.l 12008), which may show whether the 
difference between the xo of subsets is close to 90°. 

The mean and dispersion of the 28 observations allow 
us to derive a signal-to- noise map for J, Q, U and P = 
(Q 2 + U 2 ) 1 ^ 2 . For the residual systematic effects analysis, 
we have chosen a region corresponding to S/N(P)> 10 that 
includes 19 pixels. For these pixels, we compute the mean 
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Figure 2. Maps of the Crab nebula at 89.189 GHz in antenna temperature (K) for intensity I (top-left), Q polarization 
(top-right), U polarization (bottom- left) and polarized intensity P (bottom-right). The position of the Crab nebula pulsar 
is indicated by the black cross. 



polarization angle «sk y ,j = 0.5 • (arctan(C/ 7 -/Qj)) and the 
mean polarization fraction = (Pj/Ij) for each map j. 
Then, <^sky,j and are averaged for each subset to give 
the mean asky and A and the error on the mean is given by 
the standard deviation over the subset of maps, preventing 
correlations between Q and [/, or among pixels. 

Results are displayed in Table [TJ Values are stable 
against the selection of different data subsets. For the po- 
larization angle, for each subset selection criteria, the val- 
ues are compatible at less than two a, showing no clue to 
residual systematic effects in the data. We observe a sim- 
ilar behavior for the polarization fraction, where only the 
subset corresponding to a and S scans shows a departure 
of slightly more than two a. This systematic behavior is 
due to a percent-level difference in the intensity maps cor- 
responding to each type of scan, caused by residual errors 
in the destriping which is less precise for ^-scanned maps 
due to the limited number of off-source pixels in this case 
and leading to an overestimation of the intensity. 



5. Results 

After the raw data reduction, including atmospheric decon- 
tamination by linear baselines subtraction, we obtained a 
set of 28 systematic-reduced individual maps with which we 
compute average and standard deviation maps for /, Q, U , 
P = (Q 2 + /7 2 ) 1 / 2 , a Sky = 0.5 • arctan(/7/Q) and n ee P/I. 



The Crab nebula intensity / map is displayed on top- 
left panel of Figure [2j in antenna temperature, showing a 
maximum 0.91 K emission half an arcmin south from the 
pulsar position. The flux density integrated over the source 
and its error bar derived from the pixel-to-pixel standard 
deviation are measured as S{I) = 195.5±11.0 Jy. This value 
is compatible within 3a with the value given by the WMAP 
satellite at 92.9 GHz of S(I) = 229 ± 11 Jy (Pase et al. 
120071) . 

The Q and U maps of the Crab nebula are displayed on 
top-right and bottom- left panel of Figure [2j in antenna tem- 
perature (K). Q map is showing a roughly positive value, 
up to 0.16 K while U is showing a negative one down to 
-0.20 K. 

The polarized intensity P map is displayed on bottom- 
right panel of Figure [2j in antenna temperature, showing a 
maximum of 0.25 K half an arcmin south from the pulsar 
position. We can see that the Crab nebula is less extended 
in polarized intensity than in Stokes /. The polarized flux 
density is measured as S(P) = 14.5 ± 3.2 Jy, leading to a 
mean polarization fraction for the whole Crab nebula source 
ofn = 7.4±0.7%. This value is com patible with the value 
measured by WMAP of 7.6 ± 2.0 % (|Page et al.ll2007l ). 

The asky map is displayed on the left panel of Figure [3l 
It is worth noting that asky is almost constant around 150° 
in the region of maximum polarized intensity (P > 0.05 K). 
This region corresponds to the region where each pixel's 
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Figure 3. (left): Map of the polarization angle asky in degree on which contours of the polarized intensity P for 0.04, 0.08 
and 0.16 K have been overploted. (right): Map of the polarization fraction II on which has been displayed the intensity 
contour for 0.15, 0.30 and 0.45 K. 





n(%) 


^Sky (°) 


Pulsar position 


13.9 ±0.6 


158.1 


±0.5 


S/N(P)> 10 region 


23.8 ±0.5 


152.1 


±0.3 


S/N(P)> 3 region 


15.6 ±0.3 


153.7 


±0.4 


seen by a 5' beam 


8.8 ±0.2 


149.9 


±0.2 


seen by a 10' beam 


7.7 ±0.2 


148.8 


±0.2 



Table 2. Values of the polarization angle asky and of the 
polarization fraction II. 



Furthermore, we convolved each individual map by a circu- 
lar Gaussian beams of 5 arcmin and 10 arcmin of FWHM 
before computing the mean polarization angle and fraction. 
These beams mimic the Planck satellite's beams (10' for 
the 100 GHz and 5' for the 217 and 353 GHz channels) and 
give an estimation on how a generic CMB experiment hav- 
ing Planck-like characteristic would see the Crab nebula's 
polarization. For each of these cases, the mean value and 
its associated error were computed similarly to section 14.41 



standard deviation cr(asky) is the lowest (< 30°), going 
down to cr(asky) ^ 3° for the most intensely polarized pix- 
els. The value of the angle we measured at the flux peak, 
<^Sky = 149.0 ± 1.4°, is in very good agree ment with previ- 
ous measurements at 135 and 1100 /im ([Flett fc Murray! 
119911 : iGreaves et aD 12003). On the map edges, asky pixel- 
to-pixel deviation is high, except for the northern region 
where the average angle is around 230°. When looking at 
the cr(asky) map, we can see that pixels on the edges have 
an undetermined value (<r(ask y ) — 90°), whereas in the 
northern region they have 30° < cr(asky) < 50°, indicating 
a polarization angle there that is significantly different from 
that near the center of the nebula. 

The polarization fraction n map is displayed on the 
right panel of Figure [3j We have set to the pixels for 
which intensity is lower than 0.02 K in order to avoid the 
divergence of n. Maximum polarization fraction is found to 
be spatially correlated with maximum polarized intensity 
region, with a polarization fraction reaching 30% for few 
pixels. A 1.5 arcmin circular region around the maximum 
polarization fraction has a significant n ~ 20%. 

To summarize all these results, we display in Figure 
[H the polarized flux density P and the orientation of 
the polarization vectors associated with each map pixel. 
The polarization vectors pat tern is compa tible with 9 
mm (|Flett fc Henderson! 1 1979f ) and 850 //m ([Greaves et al.l 
I2003D observations. 

We finally computed several values of the polarization 
angle asky and of the polarization fraction n for differ- 
ent regions, displayed in Table [2j We selected the pulsar 
position and high signal-to-noise ratio regions for the po- 
larized intensity (S/N(P)> 3, corresponding roughly to 
the source's polarized emission extent and S/N(P)> 10). 



6. Extension to other frequencies 

The total intensity emission of the Crab nebula from 1 to 
10 6 GHz is dominated by the well known synchrotron ra- 
diation observed at radio wavelengths with only one ex- 
tra d ust component in the far infrared (|Macias-Perez et~al] 
l2Q10h . This synchrotron emission shows a decrease of flux 
with increasing frequency which can be repre sented by a 
powe r law of spectral index (3 = —0.296 ±0.06 (Ba ars et al.l 
[1977 : M acias-Pe rez et al.ll2010l ) from the radio to the sub- 
millimeter domains. The flux is als o decreasing with time a t 
a rate of a = 0.167±0.015 % yr" 1 (jAller fc Revnoldslll985f ). 
Moreover, from the visible to the X-rays the synchrotron 
emission evolves towards a much harder spectrum repre- 
sented by a power law of s pectral index j3 = —0.698 ±0.018 
(|Macias-Perez et al.ll2010l ). 

These statements allow us to postulate that where the 
synchrotron dominates, the emission at different wave- 
lengths is produced by particles accelerated by the same 
magnetic field. The direction of polarization is thus ex- 
pected to be constant while the polarization fraction may 
vary. Nevertheless, we expect values for both the polariza- 
tion angle and the polarization fraction, at low-resolution, 
to be similar in the millimeter and submillimeter where the 
emission is produced by the same electron population. 

Values of the polarization angle of the Crab nebula 
source have been reported over this wide range of wave- 
lengths from the radio (e .g. at 9 mm, asky = 154.8 ± 2.0, 
iFlett fc Henderson! 0.979) to the millimeter (at 3.3 mm, 
c^Sky = 153.7 ± 0.4, this paper) and to the X- rays (at 240 
pm, asky = 155.8 ±1.4, Weis skopf et al.| [l978). We can see 
that the value of the polarization angle is strikingly con- 
stant over nearly ten decades of wavelength. Furthermore, 
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Figure 4. Map of the Crab nebula polarized intensity P at 89.189 GHz in antenna temperature on which polarization 
vectors have been overpolted. Intensity contours at 0.15, 0.30 and 0.45 K are also displayed. 



the value measured in the X-ray, similar to the one we mea- 
sured at 90 GHz, indicates that we are probably dominated 
by synchrotron emission for the same regions in both cases 
and that an extrapolation up to 353 GHz, where present 
and future CMB experiments such as the Planck satellite 
are observing the sky, is rather safe. 

The polarization fraction can be compared to other 
measurements too, keeping in mind that it is a quantity 
which changes inside the source rapidly and that shall be 
compared only for experiments observing with a compara- 
ble beam size. The comparison of our value of n = 15.6 ± 
0.2% to radio observation s at 9 mm giving n = 16 d= 1% 
([Flett fc Henderson! 1 1979f ) and the value we obtain when 
our maps were convolved by a 10' b eam of 7.7 ± 0.2% to 
the WMAP experiment (7.6 ±2.0 %, iPage et aLlkoOfl ) in- 
dicates that those measurements are coherent with an emis- 
sion coming from the same electron population wich leads 
to a constant polarization fraction over the CMB frequency 
range. 



7. Summary 

We mapped the polarized emission of the Crab nebula using 
the IRAM 30m telescope at 89.189 GHz with an angular 
resolution of 27", with two orthogonally linearly polarized 
heterodyne receivers. The Stokes parameters were derived 
from th e auto- and cross-c orrelations using the XPOL pro- 
cedure ([Thum et al.ll2008l ). 

Observations have been undertaken during two cam- 
pains, leading to a set of 105 individual /, Q and U maps 
having changing weather conditions, sizes and centering. 
Particular care has been taken in choosing the set of indi- 
vidual maps leading to the lower level of systematics in the 
final products. Additionnaly, jack-knife tests have been car- 
ried out in order to demonstrate the robustness of the data. 



As a result, we constructed 10' x 6.7' /, Q and U coaddition 
maps of these observations centered on the Crab nebula. 

From these maps we have computed the polarized in- 
tensity P = Q 2 + U 2 , the polarization angle asky = 
0.5-(arctan(£y/(3)) and the polarization fraction n = (P/I). 
We derived our results on the polarization characteristics 
of the Crab nebula with these maps: 

— The measured flux density is S (I) = 195 ± 11 Jy and 
the polarized flux density is S(P) = 14.5 ± 3.2 Jy. 

— The polarization angle is almost constant in the region 
of maximum emission in polarization with a mean value 
of asky = 152.1 ± 0.3°. A region north to the pulsar is 
seen with an average angle of ~230°, but it does not cor- 
respond to a strongly polarized emission region. When 
seen by a 5 arcmin beam, the mean polarization angle 
of the whole source as a value of asky = 149.9 ± 0.2°, in 
good agreement with the other measurements at radio, 
millimeter and X-rays wavelengths. 

— The polarization fraction shows a maximum in a region 
south from the pulsar position, reaching 30%. When av- 
eraged by a 5 arcmin beam, the measured mean value of 
the polarization fraction isn = 8.8±0.2%. This value is 
in very good agreement with the WMAP measurements 
at 94 GHz. 

With a precision of a few tenths of degree on the polar- 
ization angle and of a few tenths of percent on the polariza- 
tion fraction, these measurements will be especially inter- 
esting in the purpose of absolute calibration of the polar- 
ization angle and the cross-polarization leakage for present 
and future polarized CMB experiments, such as the Planck 
mission, in order to accurately disentangle instrumental and 
systematic 5-modes polarization from the primordial one. 
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